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A method for mapping known cylindrical magnetovac solutions to solutions in 
torus coordinates is developed. Identification of the cylinder ends changes topology 
from R 1 x S 1 to 5 1 x 5 1 . An analytic Einstein-Maxwell solution for a toroidal mag- 
I netic field in tori is presented. The toroidal interior is matched to an asymptotically 

Q ' flat vacuum exterior, connected by an Israel boundary layer. 
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^ ^ I. INTRODUCTION 

On 
O 

Magnetic fields can be found in many astrophysical configurations. Both poloidal and 
7-h ; n 

toroidal fields can model the observed dipolar fields in magnetic white dwarfs 111. The radius 

O 



5h 



q ' of a nonrotating white dwarf can change significantly [2] when a magnetic field is frozen-in. 

The equilibrium configuration of magnetized stars with both poloidal and toroidal magnetic 
fields has been studied for some time jjj. It is known that toroidal B fields deform neutron 
stars and it has been shown this deformation acts as a source of gravitational radiation. 

Because of the general interest in astrophysical magnetic fields, this work provides a 
method for mapping known cylindrical magnetovac solutions to magnetovac solutions in 
nested tori. The cylindrical solution used as an example has an axially symmetric magnetic 
field which is mapped to tori with a toroidal magnetic field. The cylinder solution is un- 
bounded in distance from the cylinder axis. It is straight forward to map a portion of this 
solution to a set of tori. If we conclude the map with an unbounded set of tori, nothing 
further is required. But if the torus solution is matched to an exterior vacuum at a finite 
boundary, then the metric and extrinsic curvatures, as well as the magnetic field must also 
satisfy boundary conditions. The magnetic field requires a surface current density, which 
must reside in an Israel boundary layer jfj. We construct such a layer between the magnetic 
tori and the exterior vacuum. 
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The new solution has four distinct regions. The innermost region is a singular loop 
whose time history is a singular cylinder. The next region extends from the singularity 
to the vacuum boundary. This second region contains a toroidal magnetic field in nested 
tori. At the vacuum boundary there is a third region of zero thickness with a delta function 
current density supported by an Israel layer. The fourth region is an asymptotically flat 
vacuum. 

The next section discusses the map from cylinders to tori and gives details of the mag- 
netovac cylinder solution. The cylinder solution has a singular z axis. This is mapped to a 
singular cylinder in torus coordinates. Section III gives a new solution for a toroidal magnetic 
field in torus coordinates. It is shown that the Rainich conditions for an Einstein-Maxwell 
solution are satisfied. In section IV the boundary and extrinsic curvatures are presented. 
Section V treats the magnetic field junction conditions, the boundary current, and the stress- 
energy content of the Israel layer. In section VI the electromagnetic energy, Komar mass, 
and sectional curvature mass are computed. The sectional curvature mass definition has 
been extended from spheres to include tori. We show that the magnetic energy subtracts 
quasilocally from the sectional curvature mass. Section VII presents the complete metric 
and its asymptotic and inner regions. The vacuum region is asymptotically flat and can be 
compactified. We close with a Summary. Details are collected in four appendices: Appendix 
A discusses Euclidean torus coordinates. We use coordinates (r,a,8): r labels successive 

n 

tori and (at, 3) are coordinates on a single torus Appendix B analyzes geodesies in the 
plane surrounding the singular region by examining timelike geodesies paths in an effective 
potential. Appendix C covers torus null tetrad quantities. Appendix D lists electromagnetic 
field equations and the Rainich conditions. 
Conventions 

Riemann and Ricci sign conventions are 2A v .^ a p [ = A^R^p, and R a p = R u a/3u - The 
metric signature is (+,-,-,-) and the field equations are = SttT^. We use units such 
that G = c = 1. 

Magnetic field terminology 

In cylindrical or spherical coordinate systems, magnetic field B has orthogonal toroidal 
and poloidal components. A component is toroidal if it is everywhere orthogonal to merid- 
ional planes through the cylinder axis or planes through the spherical poles. If B has axial 
symmetry, toroidal field lines are circles about the cylinder axis or spherical circles of Ion- 
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gitude. A component lying totally in meridional planes is poloidal, and is orthogonal to the 
toroidal component. 

II. CYLINDER TO TORUS 

A torus is the only compact, oriented, Euler number zero, 2-dimensional surface. A 
cylinder can be made into a torus by joining its ends. 

The cylindrical 3- metric is dr 2 + r 2 d(p 2 + dz 2 , and this is mapped to torus coordinates as 
dr 2 +r 2 d(3 2 + $t 2 da 2 . (Details of torus coordinates are given in Appendix A.) Cylindrical f — > 
torus r, and cylindrical tp — > torus (3. For the cylinder ends, dz = dtda. Since — oo < z < oo, 
and < a < 2tt, the cylinder is wrapped to a torus. This identification takes R 1 x S 1 to 
S 1 x S 1 . The cylinder z axis corresponds to 3? = r + r cos (3 = 0. Because of torus topology, 
flat metric dt 2 — {dr 2 + r 2 d(3 2 + §l 2 da 2 ) is not the Minkowksi metric. 



Cylinder Solution 

The cylindrical metric we consider (Eq.22.11, m = 2, in Exact Solutions []]]) is 



g c J u l dx^dx u = r\F/b ) 2 {dt 2 - dr 2 ) - f~ 2 (F/b ) 2 d<p 2 - r\b /F) 2 dz 2 



(1) 



with F(f) = 1 + b^f 4 . This solution satisfies the Einstein-Maxwell field equations and has 
an axially symmetric magnetic field along the z axis. The Kretschmann scalar is 



R a ^ u R a 



64 



6^(3 + 12b 2 r 4 + 62b 4 r s - 108^ 12 + 63&gf" 
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Metric g cjl has a singularity along the z axis, f = 0, where the Kretschmann scalar becomes 
infinite. 

The Ricci tensor has components t, r, ip, z 



-1 



(2) 



The magnetic field is 



(3) 
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with electromagnetic invariants 



1 b A 
-F^F*" = 16 " 



F 4 



-F^F* = 0. 
2 ^ 



The metric and magnetic field can be mapped, locally, to a different metric and magnetic 
field in a frame of nested tori. After the transform, the resulting metric is g tOT given in 
Eq.©. 

III. TORUS SOLUTION 

The toroidal magnetic metric is, with a map of cylinder cordinates {f, z, ip} to torus 
coordinates {r, a, /?}, and with 3? = r + r cos j3 and F — 1 + &o3ft 4 , 



g^dx^dx" = F 2 ^dt 2 - 9£ 4 (F 2 cos 2 + F~ 2 sin 2 /5)dr 2 - F 2 5T 2 cfa 2 



(4) 



+ K 4 (F 2 - F~ 2 )r cos /3 sin /? 2drd(3 - 9£ 4 (F 2 sin 2 /5 + F" 2 cos 2 /3)r 2 ^ 2 



with volume element 

(r^ 5 F 2 ) dtAdrAdaA dp. (5) 

Note that if parameter p is introduced by F — > F/po in metric g tor , just as parameter bo 
appears in the cylinder metric, then the Ricci and Riemann tensors are multiplied by p\. 
The limit po — > takes g tor to flatness (but not to the Minkowski metric). 

The Kretschmann scalar for g tor is R a /3 flu R a,3tJ-u ~ l/3ft 12 . The singularity of g cyl is mapped 
to a singular loop, 3ft = 0, at the center of the torus solution. The spacetime history of the 
loop is a cylinder. 

Killing symmetries are static dt and axial d a . Constant time hypersurfaces for static g tov 
have unit vector 

tA = F-^%. (6) 
In torus coordinates, the Ricci tensor has components (t, r, a,/3) 



R» = 16 



WF 4 



(7) 



-1 
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Static metrics have the timelike Killing vector as an eigenvector of the Ricci tensor. For g tov 
we find 

= ( 16 wr) % 

with eigenvalue twice the magnetic energy density. 

The metric is expanded in a null basis, g l °J = l^n u + n^l u — m fJi fn u — rh ll m l , (see Appendix 
B). The only non-zero Ricci component is 

$ n = 8b 2 ary- 4 . (8) 

With Ricci tensor 

R l ™ = -2$u(Z jU n I/ + riplv + m^rh u + m^m v ) (9) 

the Ricci "square" is 

R^K = 4($ n ) V„ 
and with $ n in Eq.fjHJ) this becomes 

b 2 ' 

Equation (fTUj) shows the Ricci "square" satisfies Rainich condition (jD14b|) . 
The vacuum limit, b — > 0, is 

g^-^dx^dx" = Wdt 2 - ^ A dr 2 - ®r 2 da 2 - r 2 $l 4 dp 2 . (11) 

^tor-vac zero j^j cc j tensor and non-zero Riemann tensor. (R Q /3fiuR a/3 ^ u )tor-va,c — 192/9? 12 . 



= 32 ^i • ( 10 ) 



Electromagnetic Field 



The vector potential for the Maxwell field is 

3? 4 



A^dx" = b j- (sin/? dr + r cos (3 d(3) } (12) 



with 



K» = 4& o(yy) K^f - S ( f } Sl r) ] • (13) 

There is no electric field since F^J has no time component and there is no current since 
V V F^ T = 0. The dual Maxwell field is 

F*, v = 26 [SMS® - 5f5^] . (14) 
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The local magnetic field _B* or = F* t^ or , which satisfies V ^B^ = 0, is 

h b 2 
M K 2 F M or 5J 2 F 4 ' 

The magnetic field lines are toroidal, i.e. B^ or is everywhere orthogonal to meridional (r, 0) 

planes through the a axis. The non-zero Maxwell null tetrad component is 

( f )l =t(2V2b )^- 1 r- 2 (16) 

The relation $ ab = a 0& between the Ricci <3> a b and the Maxwell a must be satisfied, and 
such is the case for the solution given here. This is seen explicitly since 0i0i = $n given in 
Eq.fjHJ). The Maxwell invariants are 

A tor = \F^F, V = -(0 2 + 2 ) = 16^_ (17) 
jtox = l _ F ^ F * u = o. 

The Weyl tensor vanishes on a particular torus. From Eqs. (jC4|) all Weyl components are 
zero at = 1/& 2 ,, providing an effective "magnetic range". 



IV. INTERIOR-EXTERIOR METRIC AND EXTRINSIC CURVATURE 

A. Boundary 

One can match interior metric g tOT to exterior metric g tor - vac . The boundary surface is 
located at 

r = Tb = const. 

The metric match is straight forward since g tor - vac is g tor with bo = 0. The second junction 
condition requires extrinsic curvatures to match at r^. is the unit normal at r;,, N^N^ = 
-1, N^dx" = Ndr, 

N = ^ 2 r[i + b 2 Xsm 2 p(i + r)(i + r 2 )}- 1/2 , iv vac = ^ 2 . (is) 

The extrinsic curvature, K^, is the projected covariant derivative 

K% = N p , a . (19) 

where -L P u= <? p „ + NPN^ projects into the boundary. 
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B. Interior 



At r = r b , fC° r , with = r + r b cos (3 and F b = 1 + fr 2 ,^ 4 , nas components 



2JNkcos/?(l + 36g3Kg) 
* WW = ^ > (20a) 



iV fc cos/j(l-3&M 

*(a)(a) = ^ > ( 20b ) 



N b = U 2 b F b [l + b 2 X sin 2 0(1 + F b ) (1 + F fe 2 )]- 1/2 



(0)(0) ~ ~%Ff (2 ° C) 



W 6 = 3%/r 6 - 2r /r b + b 2 X(3U b /r b + 2r /r b ) - 2&X(4 cos 3 /3 + 9 cos (3 - 5r /r b ) 
- (b 6 + C)^(10^/r 6 - 28K 4 cos (3 sin 2 /? - m b /r b - 2r /r b ) 
+ 2b 8 X(5%/r b - 14 cos sin 2 (3). 

Relevant g^ T components are 

gT = ^" 4 r- 2 , gtl {a) = -^ 2 f- 2 , giT = -r- 2 5R- 4 (r 4 sin 2 p + cos 2 

The trace of K^J is 

= -m l2coaff{l + 3W) - cosi3(1 - 3 « + rnn*?e+«*0) ] - (21) 

C. Exterior 

fT*° r ~ vac has components 

K^ c = -2%cos(3, Kl^:f = -cos(3/Kl = r b M b (3M b - 2r ). (22) 

Relevant g^-vac components are 

(*)W = («)(«) _«2 (/3)(/3) _ r -2»-4 

it tor- vac JL ? if tor- vac Jl - 5 y tor- vac ' Jl - 

The trace of Kfi-™ is 

^ tor " vac = ~ (cos/? + 3%/r fe - 2r /r 5 ). (23) 
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V. MAGNETIC FIELD MATCH 

With unit normal N, the magnetic field at the boundary requires 

(B t or ■ N)b = (fivac ■ N vac ) b . (24) 

This condition is satisfied since the normal components are zero on both sides of the bound- 
ary. The tangential component must obey 

(B toI x N) b = -J s (25) 

where J s is a surface current density. The jump in tangential component (non-zero to zero) 
requires an Israel surface layer to support J s . 



Israel Layer 

If the extrinsic curvatures do not match at the boundary, an Israel boundary layer is 
created. The stress-energy content of the Israel layer is constructed from the mismatch in 
the curvatures. The stress-energy of the boundary layer is [8J 

-87rS"t= - [K] ±» v . (26) 

Here K = and [K] := i^ cxterior — /^interior d eno t es the jump of K (for any tensorial 
quantity defined on both sides of the boundary). The non-zero components of the stress- 
energy are S'S, S^, and S^}. We find 

Pg] = -g^l«1(«-/» + - W + Ft J^ ms ^ (27a) 

-iV 6 cos/3(l-36X)] 

[^)] = W3^^ 2 cos ^ + 3 ^ ~ 2r °/ r ») - rA • if" 2R (27b) 
[ ' vi-b^i Ffc sin p + cos 2 p 



-2JV 6 cos/3(l + 36gW, 



[^g] = ^ - iV 6 (l + 9&X)] (27c) 

When the magnetic field is set to zero, i.e. b$ —>■ 0, all jumps vanish. 

A graph of Sj$ for weak magnetic field, with F \ ~ 1, is bell-shaped and peaked at 
(P) = n. (Fig. 1) With large magnetic field, F \ ~ &o^&' the g ra ph has two bell-shaped 
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peaks at ((3) = 7r/2,37r/2 falling to zero at ((3) = it. (Fig. 2) The other components have 
similar graphs. The stress-energy is peaked at the "top, 7r/2" and "bottom, 37r/2" of a given 
torus. Since magnetic energy subtracts quasilocally from mass energy, the stress-energy is 
greater at top and bottom to maintain the torus shape. With no magnetic energy, the 
stress-energy is a maximum at the innermost point, (/3) = ir, to balance mass attraction. 



0.011 H 




1 2 3 4 5 6 



FIG. 1: ro/rfc = 10. Weak magnetic field with Fjf ~ 1. Vertical scale adusted to graph. 




FIG. 2: Large magnetic field. F b 2 ~ bffif. Vertical scale adjusted to graph with 6q = 1- 
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In the outer equatorial plane, with cos (3 = 1, 



~ + r /r b )(r b - 4r ) + 0(r b 5 ), (28a) 



'(t)J 

7(0)1 



~ -2^(1 + r /r 6 )(r 6 - 2r ) + 0(r b 2 ), (28b) 

'(/3)J 



[S$] ~ — 6g(l + r /r fe ) + 0(r b 5 ). (28c) 



Poisson 



expresses the complete energy-momentum tensor as 

= e(/)T- tcrior + e(-o^ teriOT + 5(i)s^. 

In this work T^ tcrior = for # tor - vac and T^ terior = T^° r for g tOT . 9 is the Heaviside step 
function and 5(1) is a Dirac distribution of shell width I. We have 

?V = 0(r - r fe )T;° r + 5(r - r 6 )<V 

The conservation equation for S^, its covariant divergence, expresses the jump in T*° r 
components at the boundary. To find the surface current density, we write the field tensor 
at the boundary 

-^interior = @( r — r &)-^tor 

where 

V^ntcrior = = K r ~ n)F^ ] . With J s from Eq.(J2SJ) the current is 

AnJ^ = 5(r-r b )(-^)d,. (29) 

Any one of the current lines is an a = const circle around the torus. Together, the circles 
form a surface current on the bounding torus. 

The graph of one of the current lines, J s vs (/?), is bell-shaped and peaked at ((3) = it. 
For weak magnetic field, with F b ~ 1, the shape is dominated by 9f^ 5 ~ (ro/r b + cos/3) -5 . 
The curve narrows and the peak rises as r b — > tq, while the volume containing the magnetic 
field shrinks, with a spike as the volume nears zero. With large magnetic field, F \ ~ bffif, 
the current goes as 9^ 13 . The curve shape is the same, with smaller magnitude. 
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FIG. 3: J s vs {(5). F 2 ~ 1, ro/rj, = 1.1 Vertical scale adjusted to graph. 



VI. MASS AND ENERGY 



A. Magnetic energy 



The electromagnetic energy follows from integrating the magnetic energy density, u, over 
t = const hypersurfaces within 3-volume v / —gd 3 x = r3? 5 F 2 dadfldr 

u = -\b^ 

= 8b 2 /(^ 2 r 4 ). 



The toroidal energy is 



r 2tt 



C/ tor (r) = 2bl J J 



r 



r'(ro + r' cos /3) 3 
[l + b 2 (r + r' cos/5) 4 ] 2 



d/3dr'. 



Termwise integration provides 

U toI (r) -- 



b 2 r 3 



1 + 



2„4\2 



(r - r ) 2 + 0(r - r f 



(30) 



B. Komar mass 



The Komar superpotential, with Killing vector /c M , is 



(31) 
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For k p = 5^ we find 



V 2 

"komar ~ U ^ 



The Komar mass is expressed by 

M(d t ) 



8tt 



(32) 



The mass integral, for i = const, r = const 2-surfaces with dS^ = 2(t^r tU — t tlJ r^)dadj3, is 



2tt 

M tOT (r) = 4r J cos/3 
o 



1 + 3^(r + rcos ( 5)" 
1 + &o( r o + rcos/3) 4 



Integrating termwise provides the toroidal quasi-local mass 



M tnv (r) 



32vr& 2 r 3 1 



r 2 + Ofr" 



The mass expansion is finite since vacuum starts at r = r&. 



d/3 



(33) 



C. Sectional Curvature mass 

Spherical fluids, with 2- metric r 2 ((ii? 2 + sin 2 $d(p 2 ), have a well-defined sectional curvature 
mass given by 

-2m sph (r) = r 3 R^ p J^ u $ p 0° 

where the Riemann tensor is weighted by r 3 because of spherical geometry. The Reissner- 
Nordstrom solution (Schwarzschild mass m , charge q) has sectional curvature mass 



2m rn (r) = 2m — q 2 /r. 



(34) 



The electrostatic field energy due to charge q has a negative quasilocal contribution. 

One can extend this notion to energy in tori. Vectors a p and f3 p span r = const tori 
and are surface-forming, i.e. P u ol^. v — a 1 ' (3 P . V = aa p + b(3 p . We define the torus sectional 
curvature mass as 

-2m tor (r) = [^R pi/pa a^ u a p ^} W)=0 (35) 

with the Riemann tensor weighted by 9ft 6 because of the torus geometry in metric (jlj). Fixing 
(/3) = selects 9ft = r + r for a given torus. We find 

l_ & 2 (ro + r) 4 



m t or(r) = m 



[i + b 2 (r + r yy 



(36) 
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Parameter m appears in the Riemann tensor when metric function f is scaled by m . 

We see that the magnetic energy has a negative contribution to m tor , just as the Reissner- 
Nordstrom electrostatic energy. In the vacuum limit, b\ ~~ > 0, mt OT (r) — > m . 



VII. COMPLETE METRIC 



Metrics g tor and g tor - vac combine to provide 



^toroidal _ J 9 tOT > r o < r < r b 
5f tor_vac , r b < r < oo 



A. Inner region 



The inner torus boundary has coordinates cos/5 = —1 and r = r , with metric functions 
3ft = 0, and F = 1. The toroidal metric goes as g tOT ~ 3ft 4 (cft 2 - dr 2 - r 2 d(5 2 ) - (l/3ft) 2 da 2 
and is not well behaved. Since BTR^ = (32b 2 /3ft 2 F 4 ) 2 and R a ^ v R a ^ v ~ l/3ft 12 , it is clear 
that g tor has a singularity at r = ro, cos/3 = — 1, the 3ft = loop. The spacetime history of 
the central singular loop is a singular cylinder. This reflects the singular axis of the original 
cylinder metric (JT|). 



B. Asymptotic region 

The vacuum region starts at r&. For large distances r >> ro, and with | cos (3 |~ 1, metric 
function 3ft = r + r cos /3 go as 3ft ~ r. 
(yf tor is matched to a vacuum solution 

g t ™-™dx tx dx v = 3ft 4 (rft 2 - dr 2 - 3ft" 6 rfa 2 - r 2 d/3 2 ). 

Asymptotically, for | cos/5 |> 

^tor-vac _ r 4^ t 2 _ ^2 _ ^6^2 _ r 2^2) (gg) 

When (/?) is near tt/2 or 37r/2 and cos/3 ~ e (near zero) 

^or-vac _ e 4 r 4^ f 2 _ ^2 _ ^6^2 _ ^2) 
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There are two directions which do not have a well defined limiting metric: (/?) = n/2, 37r/2, 



by choosing different torus coordinates. Morse and Feshbach [10] use torus coordinates with 
a hyperbolic radial coordinate. Present coordinates {r, a,/3} tor map to {//, 0,?7}mf- The 
analog of (/?) is angle (77). The Gaussian curvature of (0, 77) 2-surfaces is 



and for (77) = 7r/2, 37r/2, /Cmf = 1- 

With r 4 as a conformal factor, g tor - vac can be compactified and is asymptotically flat. In 
the compactified diagram, the boundary of the singular cylinder curves from I to I + while 
the null surfaces of the tori extend to X + . The spacetime of g tor is asymptotically flat and 
foliated by nested tori. 



We have developed a method for mapping known static cylindrical magnetovac solutions 
to static solutions in torus coordinates. The cylinder topology has been changed from R 1 x S 1 
to torus S 1 x S 1 by identifying the top and bottom cylinder ends. An explicit example has 
been presented with a toroidal magnetic field. The example satisfies Rainich conditions for 
valid electrovac solutions. 

The new solution has three parameters: r specifies the torus geometry, 6q the magnetic 
field and magnetic energy density, and mo is the sectional curvature mass parameter. The 
jump in magnetic field is supported by a surface current density in an Israel layer at the 
vacuum boundary. The solution has an asymptotically flat vacuum exterior, and an Einstein- 
Maxwell interior with a singular cylinder at the center. By extending sectional curvature 
mass to include energy in tori, we have shown that the magnetic energy subtracts quasilocally 
from m t0 r- 

There are other exact cylinder solutions, such as the Melvin metric, which have no singular 
region but whose map to torus coordinates provide a more complicated magnetic field. These 
will be studied in future work, which will also add rotation to these solutions. 
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where the Gaussian curvature of the (a, (3) 2-surface is zero. This difficulty can be removed 




/Cmf = 1 — cos 77 cosh fi 



VIII. SUMMARY 
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this work 



APPENDIX A: TORUS COORDINATES IN EUCLIDEAN 3-SPACE 

A torus can be generated by sweeping a circle, orthogonal to the x-y plane, about the z- 
axis. In the x-y plane tq is the distance from the z-axis to the center of the generating circle, 
and r is the radius of the generating circle. For ranges < r < oo, < a < 2tt, < (3 < 2tt, 
with 3? = r + r cos (3 

x = 3?cosa, y = $lsma, z = rsm(3. (Al) 
Cartesian and toroidal coordinates are related by 

2 , 2 i 2 cr>2 i 2 • 2 n 

x + y + z —M + r sin (3. 

The torus r = lies in the z = plane and has circular radius r . Looking down the z-axis 
(about which a has range < a < 2ir) at the torus, one sees the torus boundaries as two 
concentric circles. The (3 = constant surfaces, < (3 < 2tt, are spheres centered on the 
z-axis. In the plane (3 = 7r/2, x 2 + y 2 + z 2 = Tq + r 2 . 

The Euclidean metric dx 2 + dy 2 + dz 2 with coordinates (|A1|) becomes 

dr 2 + $l 2 da 2 + r 2 d(3 2 . (A2) 

The (a,/3) 2-surface has Gaussian curvature 

K = (A3) 
r(r + r cos [3) 

The angle f3 ranges from to 2n over the generating circle, while /C is negative, zero, and 
positive. 

K < for tt/2 < (3 < 3tt/2 

/C = for (3 = vr/2, 3tt/2 

/C > for < (3 < vr/2, 3tt/2 < (3 < 2ix. 

The Euler-Poincare characteristic is 

x = )C r(r + r cos /?)cf/3 = (A4) 



o 



corresponding to the torus (sphere with one handle). 
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APPENDIX B: GEODESICS 

To better understand the region surrounding the central singularity we examine the 
geodesic flow. From metric g tOT ~ vac we have a Lagrangean for geodesic paths 

C = h^H 2 - K 4 f 2 - 5T 2 d 2 - r 2 ^ 4 /3 2 ) (Bl) 

where overdots denote d/ds. The geodesic equations follow from 

d_oc _oc =Q (m) 

dsdx a dx a K ' 

Killing symmetries d t and d a yield first integrals 

—j = E = 5ft 4 t, 
dt 

d£ t <rv-2 • 

— = -Jo = -3? a 
oa 

or 

i = £ /^ 4 , a = J 3? 2 (B3) 
For x a = ((3), the timelike geodesic equation is 

-r 2 $l 4 f3 + 2r 2 ^ 3 (rsmp)p 2 + 2K 3 (r sin (3){i 2 - r 2 ) + 3? 3 (r s\nf3)a 2 = 0. (B4) 

If {(3) = 7i then j3 = 0. Furthermore, if /3 — and /3 = 0, then the orbit remains in the 
((3) = 7i plane. (In order to examine the singular region, we choose ((3) = n rather than 
{(3) = 0.) Again, from metric g tor - vac 

1 = E 2 /^ - J 2 K 2 - K 4 r 2 - r 2 $l 4 $ 2 

1 = E 2 /(r - r) 4 - J 2 (r - r) 2 - (r - r) 4 r 2 (B5) 

We follow the central force problem of classical mechanics, where E = \mr 2 + V e g(r), 
and write r 2 + Kfr( r ) = const. Thus, the timelike geodesic paths in the {(3) = n plane can 
be described by the effective potential 

V eS (r) = l/(r - r ) 4 + J 2 /(r - r ) 2 - E 2 /(r - r ) 8 . (B6) 

For Jq = E 2 , the graph of V cS has a hard negative core at r = r . The graph increases 
exponentially to positive values with a maximum at r max , and then falls to zero as 1/r 2 . 
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FIG. 4: V eS with r = 0, J 2 = E% = 1 

r max is a positive root of Jq (r — r ) 6 + 2(r — r ) 4 — 4Eq = 0. At r = r max there is a closed, 
unstable, geodesic path around the singular region. 

For Jq = 0, the shape of the graph is the same, but the height of r max is decreased. 

r max = r + (2£ 2 ) 1/4 - 

For large angular momentum, Jq > Eq, the graph has no negative values. 

APPENDIX C: TORUS NULL TETRAD QUANTITIES 

The metric, g*° r = l^n u + n^l u — m^m u — rh^m u with 3? = r + r cos (3 and F — 1 + fr 2 ^ 4 , 
has the null basis 



l^dx" = ■^(Vfdt-M^da), 

n^da? = -^=($l 2 dt + ft- 1 da), 
V2 

3? 2 

m^dx^ = i=[{F cos (3 + iF 1 smP)dr — (F sin/3 — %F 1 cos 0)rd(3], 

V2 

K 2 

m^dx^ = i=[{F cos 13 — if 1 sin (3)dr — (F sin f3 + iF 1 cos f3)rdf3}. 

v2 



(Cla) 
(Clb) 

(Clc) 

(Cld) 
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The spin coefficients are 



e = a = A = p = // = 7 = (C2a) 

K = ~^ (C2b) 

i/ = -k (C2c) 

7r = J_f_J_) (1 + 96X) ( C2d ) 



2V^ 3^ 3 r 2 

r = -7T (C2e) 

« = -^(^)(!- 9 ^ 4 ) (C2f) 

/3 = -a (C2g) 

m M and m M are null geodesies, since m^ v m v — {(5 — a)m M . The only non-zero Ricci 
component is 

*"=«^ «*> 

The Weyl tensor null tetrad components are 

*o = -(|^)(l -&§»*) (C4a) 

= ^ 3 = (C4b) 

*2 = -(-^)(l-^)(l + ^) (C4c) 

^ 4 = (C4d) 

which correspond to Petrov type I. 

The quadratic Ricci and Weyl invariants are 

W V R^ = (326 2 /K 2 F 4 ) 2 , (C5a) 

rafi^r 192 (l-6g^)(l + 66g^ + 216ggy) 

W^- K12/r4 (1 + 66^ + 6^8) • 

APPENDIX D: RELATIVISTIC ELECTROMAGNETIC FIELD 

The Maxwell field is represented by skew tensor F^ U) with dual field 
F»v = \\f~9 e ' iivafi F al3 . The symmetric energy-momentum is 

T, u = ^{F m FZ + \g„ v F a pF^). (D6) 
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In a local Lorentz frame, the energy density is 

T 00 = -(E 2 + B 2 )/8n. (D7) 

The Einstein field equations provide the trace-free Ricci tensor 

R P u = ~{F, a F a v + ^FopF**). (D8) 

For static systems there exists a hypersurface orthogonal timelike Killing vector = d t . 
The constant time hypersurfaces form a family of Lorentz frames with unit normal i v = 
£ u / i^n^) 1 ^ 2 ■ Local electric and magnetic fields are defined by 

Ffx = F,J» = (0, E), B,= F*J V = (0, -B). (D9) 

Maxwell's static equations are 

V • E = Anp, Vx£ = 0, V • B = 0, V x £? = An J. (D10) 

The field invariants are 



h := \f„ v F^ = B 2 - E 2 , I 2 := ^F^F^ = -2E^. (Dll) 



The Maxwell components on an anti-self dual bivector basis are 

~(F» V + iF*^) = <f> U f " / - <piM^ + (p 2 V^. (D12) 
2 

The field invariants are 

h +ih = 4(0o02 - </>?)• 

There are three necessary and sufficient Rainich conditions for to be an electrovac 
energy-momentum tensor. 
1. T^y (and Ricci) must be trace-free. 

g^T^ = g^R pv = 0. (D13) 



2. The "square" of Tn V is a positive multiple of the unit tensor. 



= (A 2 + ll)K (D14a) 
or R^R% = ^(R a/3 R a ^ u . (D14b) 
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3. Since the electromagnetic energy- momentum density, T 00 , is negative in a local Minkowski 
frame, it is necessary for any timelike vector that 

T^tH v < 0. (D15) 
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